Brownian vortexes 
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A particle diffusing around a point of stable mechanical equilibrium in a static but non- 
conservative force field enters into a steady state characterized by circulation in the probability 
flux. Circulation in such a Brownian vortex is not simply a deterministic response to the solenoidal 
component of the force, but rather reflects an interplay between force-driven probability currents 
and diffusion. As an example of this previously unrecognized class of stochastic heat engines, we 
consider a colloidal sphere stably trapped in a conventional optical tweezer. Rather than coming 
into thermodynamic equilibrium with the surrounding heat bath, the particle's Brownian fluctua- 
tions are biased into a toroidal roll. We demonstrate both theoretically and experimentally that the 
circulation in this practical realization of the Brownian vortex can undergo flux reversal. 



Stochastic heat engines such as thermal ratchets and 
Brownian motors use non-conservative forces to eke a 
flux of energy or probability out of otherwise random 
thermal fluctuations [lj]. Unlike deterministic machines, 
these systems would not operate at all without the mo- 
tivating influence of thermal noise. Although a great 
many mechanisms have been proposed for implement- 
ing stochastic heat engines, virtually all rely on time de- 
pendence in the force field to rectify fluctuations. Here, 
we demonstrate that time-independent force fields also 
can create stochastic heat engines, provided the force 
has a non- vanishing solenoidal component. The result- 
ing steady-state probability currents feature circulating 
rolls, which inspire us to name such systems Brownian 
vortexes. We demonstrate both theoretically and exper- 
imentally that Brownian vortexes arise naturally in the 
motions of colloidal spheres localized in optical traps, and 
that the probability flux can reverse direction with con- 
tinuous changes in such control parameters as the laser 
power or temperature. 

Our discussion focuses on the behavior of a single par- 
ticle immersed in a viscous medium that also acts as a 
thermodynamic heat bath. The particle is acted upon by 
a static force field, 



F(r) = -V77+ V x A, 



(1) 



that is factored uniquely into a conservative irrotational 
component described by the scalar potential U (r) and a 
non-conservative solenoidal component derived from the 
vector potential A(r). In the absence of thermal forcing, 
the particle is assumed to come to a stable mechanical 
equilibrium at a point rg for which jF(rrj) = 0. 

Were it not for thermal forces, the particle would re- 
main motionless at Tq. Thermal fluctuations enable the 
particle to explore the force field with a probability den- 
sity p(r) that does not depend on time in steady state. 
The local probability flux, 



j(r) = n p(r) F(r) - y, k B T Vp(r), 



(2) 
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FIG. 1: (Color online) Schematic representation of a colloidal 
sphere undergoing toroidal circulation in an optical tweezer, 
(a) The Rayleigh limit, a < A. (b) The Mie and ray-optics 
limits, a > A. 



ature T. Because the probability is conserved, V • j = 0, 
and the flux either vanishes in equilibrium (V x A = 0) 
or else forms closed loops in steady state (V x A / 0). 
Circulation around these loops is measured by 



V x j = fipuj + [xWp x F, 



(3) 



reflects both the particle's response to the force field 
through its mobility, /i, and also its diffusivity at temper- 



where u>(r) = — V 2 A is the vorticity in F(r) ; 2]. Vor- 
texes can form in j{r) even if F(r) does not have closed 
loops because the diffusive flux represented by the second 
term in Eq. ([3]) provides the return flow. 

Both the magnitude and the direction of w(r) can vary 
with position. The rate and direction of the particle's 
circulation therefore depend on the domain over which 
the particle can diffuse at temperature T. Changing the 
temperature changes this range and therefore can reverse 
the sense of the overall circulation. The possibility of 
temperature-dependent flux reversal distinguishes Brow- 
nian vortex circulation from the more familiar interplay 
of advection and diffusion in such systems as the electric 
current flowing through a battery-powered circuit. 

As a concrete example of the Brownian vortex phe- 
nomenon, we consider the motions of a colloidal sphere 
in a optical tweezer Q, a single-beam optical gradient 
force trap that consists of a strongly focused beam of 
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FIG. 2: (Color online) (a) Streamlines of the optical force 
F(r) in the (r, z) plane exerted on a 0.4 /im radius polystyrene 
sphere in water by a 532 nm optical tweezer focused with a 
numerical aperture 1.4 lens, (b) Streamlines of the irrota- 
tional component of the force, — VU . (c) Streamlines of the 
solenoidal component of the force, V x A. Shading indicates 
the relative probability density p(r) at T — 21°C and a laser 
power of 0.3 W. 

light [|. Recent particle-tracking measurements [4j over- 
turned the conventional understanding that a trapped 
particle comes to thermodynamic equilibrium within an 
optical tweezer, revealing instead a circulating steady 
state indicated schematically in Fig. [TJa). The original 
explanation for these observations |4J relied on an ideal- 
ized model for the optical forces acting on a sphere in 
which intensity gradients drew the sphere to the focal 
point with a Hookean force and radiation pressure dis- 
placed the particle downstream with a force proportional 
to the local intensity. This model predicts the appearance 
of a single toroidal vortex in p(r) directed downstream 
along the optical axis and returning upstream at larger 
displacements. 

Light scattering by colloidal particles is more compli- 
cated than suggested by this model, particularly for par- 
ticles such as those studied in Ref. [4[ that are larger 
than the wavelength of light [!, 0] . Figure EJa) shows 
streamlines of the force field computed with Lorenz-Mie 
theory [3] for a micrometer-diameter polystyrene sphere 
in water trapped within an optical tweezer. The trap is 
formed from a laser beam with a vacuum wavelength of 
A = 532 nm propagating in the z direction and brought 
to a focus by an ideal lens with numerical aperture 1.4. 
Streamlines are projected into the (r, z) plane in cylin- 
drical coordinates, r — (r,(f>,z). Figures [^b) and (c) 
show the irrotational and solenoidal components of the 
force field, respectively, and were obtained through the 
Hclmholtz- Hodge decomposition [2| of F(r). 

The particle's thermally driven trajectory r p (t) 
through this force field was computed with a Brownian 
dynamics simulation of the Langevin equation 

r p = pF(r p )+pf(t), (4) 

where // = (67r^a) _1 is the Stokes mobility for a sphere of 
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FIG. 3: (Color online) Computed streamlines of a 0.4 /urn 
radius polystyrene sphere diffusing in the optical force field of 
Fig- [2] (a) Forward circulation at P = 0.7 W. (b) Concentric 
counter-rotating vortexes at P = 0.5 W. (c) Complete flux 
reversal at P = 0.3 W. 

radius a moving through a fluid of viscosity 77, and where 
the stochastic force due to thermal fluctuations satisfies 
(f(t)) = and (f(t) ■ f^)) = 2[i- 1 k B T6(t - t'). The 
probability distribution and flux then are computed with 
nonparametric density estimators Q as 

*■> = (^(-"(-^9^)) » d < 5 > 

where a is chosen to minimize the variance in p(r) with- 
out unnecessary blurring. This approach lends itself to 
direct comparison between simulation and experimental 
measurements of colloidal trajectories. 

Shading in Fig. [5] represents the probability density 
p(r) at T = 21°C for an optical tweezer powered by 
P = 0.3 W. Neither, F(r) nor its solenoidal component 
display loops in this range. Nevertheless, the streamlines 
°f j( r )i plotted in Fig. [3] show loops in the (r, z) plane 
consistent with the appearance of a toroidal vortex cen- 
tered on the optical axis. It should be emphasized that 
the observed circulation is apparent only in a very long 
trajectory or in an ensemble average of shorter trajec- 
tories; at short times, the particle's motion resembles a 
random walk in a harmonic well, its toroidal bias too 
subtle to be perceived. 

Figure [3fa) shows streamlines at a comparatively large 
laser power, P = 0.7 W, for which the particle is well 
localized near the optical axis. Under these conditions, 
the particle circulates in a single toroidal vortex, much as 
was predicted in Ref. [4J and portrayed in Fig. [IJa). The 
local circulation rate, [V x j(r )] ■ <j>, is uniformly positive. 

Reducing the laser power does not change the structure 
of the force field, but reduces its overall magnitude. This 
is equivalent, therefore, to increasing the effective tem- 
perature. Doing so increases the range over which the 
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particle can wander and enables it to populate a second, 
concentric counter-rotating vortex, as plotted in Fig.[3jb) 
for P = 0.5 W, and indicated schematically in Fig. QJb). 

At still lower laser power (or higher temperature), the 
outer vortex supplants the inner vortex, and the prob- 
ability current circulates once again in a single toroidal 
roll, but with its direction reversed. This flux reversal is 
demonstrated in Fig. life) for P = 0.3 W. 

The origin of the flux reversal in Brownian vortex cir- 
culation can be seen in the vorticity of the force field in 
Fig. El Although the solenoidal component of the opti- 
cal force, V x A, is directed uniformly upward, V x F 
changes direction with distance from the optical axis. So 
long as the particle's probability density is concentrated 
in regions where u)(r) ■ <fi is positive, as is the case in 
Fig. Efa), the overall circulation of the probability flux 
also is positive. When the particle wanders into regions 
of negative vorticity, it circulates in the retrograde di- 
rection, as shown in Fig. [Ub). In both cases, the non- 
conservative part of the optical force field redistributes 
p{r) downstream of the beam's focal point and diffusion 
provides the return current. The single-roll structure re- 
asserts itself in the flux-reversed state when gradients in 
p(r) become large enough for diffusion to outstrip advec- 
tion along the optical axis. 

We also observed flux reversal in Brownian vortex cir- 
culation through experimental observations of colloidal 
spheres trapped in optical tweezers. Our system con- 
sists of 1.5 nm diameter colloidal silica spheres (Bangs 
Labs, Lot SS04N/5252) dispersed in a 50 ^m thick layer 
of water that is hermetically sealed between a clean glass 
slide and a No. 1.5 cover slip. The sample is mounted 
on the stage of an inverted light microscope (Nikon TE 
2000U) where it is observed with a 100 x numerical aper- 
ture 1.4 oil immersion objective lens (Nikon Plan Apo). 
The same objective lens is used to focus four holographic 
optical tweezers [H, [l(J EH EH arranged at the corner 
of a square with 30 /jm sides near the midplane of the 
sample. These traps are powered by a single laser (Co- 
herent Verdi 5W) operating at a vacuum wavelength of 
A = 532 nm that is imprinted with a computer-generated 
hologram [ll[ by a liquid-crystal spatial light modulator 
(Hamamatsu X7665-16) before being projected into the 
sample. The trap array is designed so that two of the 
traps have nearly the same intensity, the third is slightly 
brighter, and the fourth brighter still. This enables us 
to seek out intensity-dependent differences in simultane- 
ously acquired data sets, and thus to avoid possible sys- 
tematic effects due to vibrations or other instrumental 
fluctuations. 

The trapped spheres' three-dimensional motions are 
measured with nanometer resolution through quantita- 
tive analysis [IH of i mag es obtained with holographic 
video microscopy [3, 15 1. Holographic images are ob- 




tained by illuminating the sample with the collimated 
beam from a HeNe laser (Uniphase 10 mW) operating at 



FIG. 4: (Color online) Streamlines of the measured circu- 
lation of colloidal silica spheres trapped in optical tweezers, 
(a) k± = 6.56 pN//im: Single roll circulating in the positive 
sense, (b) k± = 2.43 pN//im: Concentric counter-rotating 
rolls, (c) k± = 2.27 pN/fim: Flux reversal. Single retro- 
grade roll. Background images show the particles' measured 
probability densities. 



632.8 nm. Light scattered by the particles interferes with 
the unscattered portion of the beam in the microscope's 
focal plane to create an in-line hologram that is magni- 
fied and recorded by a video camera (NEC TI-324AII) 
at 30 frames per second. Radiation pressure due to the 
nW//im 2 intensity of the imaging beam is negligible com- 
pared with thermal forces and forces due to the optical 
trap and so does not affect the particles' trajectories. 

Both fields of each interlaced holographic video frame 
were analyzed with Lorenz-Mie light scattering theory 
[B|, EH to measure each sphere's three-dimensional posi- 
tion r(t), with 3 nm in-plane resolution and 10 nm ax- 
ial resolution [13j A total of thirty two-minute-long tra- 
jectories were acquired at constant laser power for the 
four particles. At the end of each acquisition period, 
the trapped particles were moved automatically out of 
the field of view to acquire background holograms and 
to confirm the system's stability. Each particle's trajec- 
tories were analyzed with Eqs. |5]) and © to visualize 
the mean circulation and the results combined into maps 
of the mean circulation for each trap. In all, more than 
100,000 holograms were analyzed for each trap. 

Figure^] shows streamlines of the trajectories for three 
of the four particles, the fourth serving as a control for 
Fig. SJc). These results confirm not only the presence of 
toroidal circulation in the particles' motions, but also the 
appearance of flux reversal as a function of trap strength. 
Each trap-particle combination is characterized by its ap- 
parent 4] in-plane stiffness, k±, which is obtained from 
statistical analysis of the particle's measured in-plane 
fluctuations [11(, stiffer traps generally corresponding to 
higher laser power and lower effective temperature. The 
traps' apparent axial stiffness is a factor of 5 smaller than 
their lateral stiffness because axial intensity gradients are 
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correspondingly weaker 0, [IB] ■ 

The stiffest trap, shown in Fig. BJa), concentrates its 
particle's probability density p(r) closest to the optical 
axis and displays a single roll circulating in the positive 
sense. The weaker trap in Fig. [3(b) allows the trapped 
particle to wander further afield, where it enters into con- 
centric counter-rotating rolls, similar to the simulated 
results in Fig. [3(b). The weakest traps, one of which is 
represented in Fig.[4jc), both display a single retrograde 
roll in j(r), and thus demonstrate complete flux reversal. 

The probability distribution is centered lower in the 
weaker traps because of gravity acting on the silica 
spheres, whose 1.9 g/cm 3 density exceeds that of the sur- 
rounding water. This additional conservative force does 
not directly contribute to the particles' circulation but 
does affect what region of the optical force field the par- 
ticle occupies for a given laser power. Undoubtedly, this 
influenced the trend in Fig. [4j but does not change our 
interpretation of the phenomenon as flux reversal in a 
Brownian vortex. 

Instrumental fluctuations cannot account for our ob- 
servations because all four measurements were performed 
simultaneously in a static array of optical traps derived 
from the same laser beam. The particles are sufficiently 
separated from each other and from the walls of their 
container that hydrodynamic coupling also is unlikely to 
have influenced their motion 17|. Rather, Brownian vor- 
tex circulation, including power-dependent flux reversal, 
appears to be an inherent aspect of the statistics of col- 
loidal spheres in optical tweezers. 

In this study, we have introduced the Brownian vor- 
tex as a distinct class of noise-driven machines. Unlike 
stochastic heat engines driven by time-dependent forces, 
Brownian vortexes are driven out of equilibrium by the 
solenoidal component of a static force field. Because one- 
dimensional force fields have no solenoidal component, 
the Brownian vortex has no one-dimensional manifesta- 
tion. Not any static force field, furthermore, can support 
a Brownian vortex. For example, a force field lacking 
a stable equilibrium point cannot establish the requi- 
site probability-conserving steady-state. Still other force 
fields establish circulating steady states without thermal 
noise. An example of this is provided by the ring-like 
optical trap known as an o ptic al vortex [1 81 ] that exerts 
torques on trapped objects [191 ] through its helical wave- 
front structure [20]. These are deterministic machines 
rather than stochastic heat engines, and so are not Brow- 
nian vortexes. 

Although the illustrative simulations and experiments 
presented here focus on colloidal circulation in optical 
tweezers, the Brownian vortex is a more general phe- 
nomenon. Seeking its signature in such contexts as bio- 
logical networks and financial systems, as well as in new 
mechanical models, should provide opportunities for fu- 
ture research. Further work also is required to elucidate 
Brownian vortexes' thermodynamic properties, particu- 



larly the considerations that determine their thermody- 
namic efficiency. 
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